Classification and Stability of Phases of the Multicomponent One-Dimensional 
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The classification of the ground-state phases of complex one-dimensional electronic systems is con- 
sidered in the context of a fixed-point strategy. Examples are multichain Hubbard models, the 
Kondo-Heisenberg model, and the one-dimensional electron gas in an active environment. It is 
shown that, in order to characterize the low-energy physics, it is necessary to analyze the perturba- 
tive stability of the possible fixed points, to identify all discrete broken symmetries, and to specify 
the quantum numbers and elementary wave vectors of the gapless excitations. Many previously- 
proposed exotic phases of multichain Hubbard models are shown to be unstable because of the 
"spin-gap proximity effect." A useful tool in this analysis is a new generalization of Luttinger's the- 
orem, which shows that there is a gapless even-charge mode in any incommensurate N-component 
system. 



I. INTRODUCTION 

The basic theory of the low energy physics of the 
interacting one dimensional electron gas (1DEG), both 
with and without spin, has been well established for 
two decades. The purpose of this paper is to extend 
this general analysis to obtain a classification of the sta- 
ble fixed points of multicomponent one-dimensional elec- 
tronic systems. Examples of such problems include one- 
dimensional metals with several bands crossing the Fermi 
surface, such as multichain Hubbard ladders, |l]-^| and 
the "1DEG in an active environment," |J of which the 
most studied example is the Kondo-Heisenberg model, 
HQ i.e. a 1DEG interacting with a periodic array of 
localized spins. While these models are still one dimen- 
sional, and are amenable to the same methods of solution 
as the 1DEG, their added richness brings in significant 
new physics. In particular, in the context of the theory of 
high temperature superconductivity, this class of models 
includes some in which a spin gap and a strongly diver- 
gent superconducting susceptibility derive from purely 
repulsive interactions. Moreover, in these cases, the driv- 
ing force for the superconductivity is a lowering of the 
kinetic energy. |||| 

In ID, even at zero temperature, states with a bro- 
ken continuous symmetry are destabilized by quantum 
fluctuations. However, there are states with quasi-long- 
range order which can be characterized by the existence 
of "quasi Goldstone modes" , i. e. gapless collective modes 
of the system with a sound-like spectrum. The canoni- 
cal example of a quasi Goldstone mode is the longitudi- 
nal sound mode of a harmonic chain, For the simple 
1DEG, the relevant continuous symmetries are spin ro- 
tation invariance (SU(2)), global gauge invariance (U(l)) 
associated with charge conservation, and Galilean invari- 
ance. The latter is not an exact symmetry for the electron 



gas on a lattice but, so long as the electron density is in- 
commensurate with the crystal, the low energy dynamics 
possess an exact translational (chiral) symmetry. 



A. Classification of phases 

It is a remarkable feature of ID that all the prop- 
erties of such systems, including fermionic correla- 
tion functions, can be expressed in terms of bosonic 
fields (bosonization) corresponding to the quasi Gold- 
stone modes. Thus, it is possible to classify all 
thermodynamically-distinct ground-state phases of any 
multicomponent 1DEG by specifying: (1) any sponta- 
neously broken discrete symmetries, such as the lattice 
translation symmetry or parity, and (2) the number and 
quantum numbers of the fundamental gapless modes. 
The minimal quantum numbers of the gapless modes are 
charge, spin, and (crystal) momentum. Our convention 
will be to focus on spin and charge modes with the small- 
est non-zero momentum. Here "spin modes" have spin-1 
and charge 0, and "charge modes" have spin and charge 
2me, where m is an integer. Simply counting gapless 
modes is insufficient; for instance, a state with one gap- 
less charge and one spin mode with the same momentum 
(which we label [c, s;2kp] or [cs] for short) is distinct 
from the state in which they have different momenta, 



[c- ) 2k { p ) ][s]2k K p>] (or [c][s] for short). 

This scheme differs from the traditional method fl(J 
of classifying phases of the 1DEG in terms of the most 
divergent susceptibilities, which is appropriate when the 
goal is to understand the properties of quasi-one dimen- 
sional systems, since in most cases, weak higher dimen- 
sional couplings will stabilize a true broken symmetry 
state with the corresponding order at finite temperature. 
However, in the context of the 1DEG per se, the Lut- 
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tinger exponents, and hence the exponents governing the 
divergence as T — > of the various susceptibilities, vary 
continuously with parameters. For example, one can pass 
from a region in which the superconducting susceptibility 
is the most divergent to one in which the CDW suscep- 
tibility is the most divergent without encountering any 
thermodynamic singularities. Thus, in the strictly one- 
dimensional context, the present classification is more ap- 
propriate and, in this respect, it extends and corrects the 
ground-breaking work of Lin, Balents and Fisher jll] on 
this subject, while expanding on our earlier work on the 
general problem of a "1DEG in an active environment." 



B. Fixed-point strategy 

The concept of a fixed point of the renormalization 
group equations of a field-theory Hamiltonian was intro- 
duced by Wilson for the study of critical phenom- 
ena. This idea ultimately made its way into many-body 
theory, where the renormalization group had been used 
for some time. The main point is that the low-energy, 
long-distance physics of a given model is controlled by 
the properties of the relevant stable fixed point or crit- 
ical fixed point of the renormalization group flows. A 
particularly effective way of determining this behavior 
is to identify an exactly-solvable field-theory model that 
starts in the neighborhood of an unstable fixed point and 
flows to the same fixed point as the model in question. 
This strategy justifies the use of the "Toulouse limit" to 
solve the single-channel Kondo problem, jl3| and of field- 
theory solutions of the 1DEG with attractive backward 
scattering or umklapp scattering. jl0| It is important to 
note that a field-theory model that does not exhibit spin- 
rotation invariance may flow to a spin-rotation invariant 
fixed point Hamiltonian. The flexibility in the choice of 
solvable models allowed by this behavior is frequently ex- 
ploited, and it will be used later in this paper. 

A fixed-point strategy is more difficult to implement 
when there are many degrees of freedom, as in multicom- 
ponent systems. There may be several stable or critical 
fixed points, and it is necessary to carry out a "global 
renormalization" Jl4| in order to determine which one 
controls the low-energy physics of a given model. Usu- 
ally, such a procedure must be carried out numerically. 
Also, it is necessary to do a different stability analysis 
for each fixed point separately. The scaling dimension of 
any given operator is generally different at different fixed 
points, so an interaction that is irrelevant at one fixed 
point may become relevant at another. 

One of our major findings is that "the spin-gap prox- 
imity effect" qj is a powerful force for destabiliizing many 
putative fixed points of multicomponent one-dimensional 
systems, and enhancing superconducting correlations. 
The physics, which is driven by singlet pair tunnelling, 
is analogous to the proximity effect in superconductiv- 



ity. It serves to lock the superconducting phases of two 
subsystems, and either generates a spin gap in both sub- 
systems, or transfers a spin gap from one subsystem to 
another. We have proposed this effect as a mechanism 
of high temperature superconductivity. A more detailed 
discussion is given in Sec. III. 



C. Outline of Paper 

The paper is organized as follows: In Sec. II the the- 
ory of the IDEG is reviewed from a statistical mechanical 
point of view, and the perspective required for the anal- 
ysis of the general problem is developed. In Sec. Ill 
it is shown how this analysis can be extended to mul- 
ticomponent systems. Specifically, the classification of 
the possible fixed points in terms of their spectrum and 
broken symmetries is discussed, and the conditions for 
perturbative stability are derived. This section also con- 
tains a discussion of generalizations of Luttinger's the- 
orem to ID, including a review of a recent proof |l5| l 
of the existence of a gapless neutral collective excitation 
with momentum 2k* F , and a new theorem concerning the 
existence of a gapless charge 2me mode, where m is an 
appropriate integer. The classification of phases and the 
proof of the new theorem makes use of the continuum 
representation of the so-called 77-pairing operator, which 
is a product of fcrmion creation operators with the same 
momentum. ]l6| ] 

Readers who are familiar with the theory of the IDEG 
and are primarily interested in the illustrative examples 
may prefer to proceed directly to Sec. IV, where this 
general scheme is applied to the analysis of the global 
phase diagram of the multichain Hubbard model and 
the Kondo-Heisenberg model. Specifically, it is shown 
that I) many of the proposed partially-gapped phases 
p"T[ of the multichain Hubbard ladder are destabilized 
by the "spin-gap proximity effect" ; Q 2) the charge or- 
dered "stripe" structures which have been observed |J] 
in numerical studies of multichain t — J models have a 
(possibly non trivial) relation to the value of 2k* F de- 
rived from the generalized Luttinger's theorem ]l5{ |; and 
3) there are several, thermodynamically-distinct spin-gap 
phases of the Kondo-Heisenberg model. |L7l] 



II. THE ONE COMPONENT IDEG 

To begin with, we consider the (well understood) the- 
ory of the single component IDEG with spin from the 
perspective of its quasi Goldstone modes. All known 
zero-tempcrature thermodynamic states of the IDEG can 
be identified by: (1) spontaneously-broken discrete sym- 
metries, such as parity, (2) whether or not there is a gap- 
less spin and/or charge collective mode, (3) the smallest 
non-zero wave vector at which these modes are gapless, 
which in analogy with Fermi liquid theory is called 2k p, 



2 



(Equivalently, 2k f characterizes the long-distance oscilla- 
tory behavior of appropriate correlation functions.) (4) 
the velocity, u, and "Luttinger exponent," K, of each 
gapless mode. 



A. Definition of Terms 

The definition of the terms used above requires some 
discussion since, in much previous work, the identifica- 
tion of modes is derived from a particular calculational 
scheme rather than from general principles. Because of 
the absence of spontaneously broken spin-rotational sym- 
metry, it is possible even in the thermodynamic limit to 
classify all states by their spin quantum numbers. Thus, 
by a "spin mode" , we mean a neutral excitation with 
spin 1 ; the existence of a gapless spin excitation will typ- 
ically show up as an asymptotic power-law behavior of 
the spin-spin correlation function, (S2k F (t) ■ S-2k F (0)) ~ 
|£|i-2(5s_ (Here, 6s, is the scaling dimension of this oper- 
ator.) A charge excitation refers to excitations made by 
adding a small number of electrons to the system. Typi- 
cally, this means |lS[] an excitation with spin-0 and charge 
2e produced by adding a singlet pair of electrons with to- 
tal momentum 2kp, which will show up in the asymptotic 
form of the "^-pairing" operator, |L6| defined below. It 
is implicit in the above classification scheme that a state 
without a gapless spin mode has a "spin gap" , and that a 
state without a gapless charge mode has a "charge gap" . 

The physics of 2kp is central to the following consid- 
erations. A remarkable theorem of Yamanaka et al, |l5| ] 
which we refer to as the "generalized Luttinger's theo- 
rem," fixes 2kF — 2~Knjg at the same value as it would 
have in a non-interacting electron gas; here n is the total 
electron density and g = 2 because electrons have spin 
1/2. The theorem, which we will discuss in more de- 
tail in Sec. II, implies the existence of a gapless neutral 
excitation with momentum 2kp as long as the 1DEG is 
not commensurately locked to the underlying crystalline 
lattice. We also show that the ^-pairing operator, with 
the exact same value of 2kp, creates a gapless excitation 
in the same circumstances. When there are both gap- 
less spin and charge modes, fcj? can also be identified as 
the location of a non-analyticity in the single-particle oc- 
cupation probability, < n& >, but when there is a spin 
or charge gap, there is no sharp structure in the single- 
particle spectral function at all. Thus, one must be care- 
ful in thinking of fc^ as a Fermi momentum. 

Throughout this paper we will distinguish between 
the "excitations" and the "modes" of the system. By 
a gapless "excitation" with given quantum numbers we 
mean a set of excited states with energies which approach 
that of the ground state in the thermodynamic limit. A 
"mode" refers to an elementary excitation, which there- 
fore has a well defined dispersion relation. For example, 
above, we talked about the spin-1 charge-0 excitations of 
the system, although, in one dimension, under a broad 



range of circumstances, the elementary excitations are, 
in fact, "spinons" with spin-1/2 and charge 0. What this 
means is that the spin dynamic structure factor will ex- 
hibit a branch cut, corresponding to a continuum of two 
spinon excitations, rather than a pole corresponding to a 
magnon mode. When classifying states in terms of gap- 
less excitations, we have chosen not to distinguish which 
are elementary and which are multi-particle excitations. 

With this distinction in mind, the definition of the col- 
lective mode velocities is obvious. The existence of Lut- 
tinger exponents, K s and K c , is one of the triumphs of 
bosonization; they dictate the relation between correla- 
tion functions expressed in terms of the original electronic 
variables, and eigenmodes of the underlying bosonic free- 
field theory. For any given lattice model, the quantum 
critical exponents K a must be determined by carrying 
out a global renormalization JTifl to the appropriate fixed 
point and matching to the continuum theory In gen- 
eral, this procedure must be implemented numerically, 
either by studying the long-distance behavior of correla- 
tion functions, or by studying the finite-size scaling be- 
havior of energy levels. 

Two of the physically most important low energy fluc- 
tuations of the 1DEG are the 2kp CDW fluctuations 
and the zero-momentum (BCS-likc) pairing fluctuations. 
CDW fluctuations are neutral and spinless. For exam- 
ple, for the repulsive U Hubbard model with a half-filled 
band, the 2kp density-density correlation function has 
the same asymptotic (power-law) form as the 2kp spin- 
spin correlation function, although there is manifestly a 
charge gap. p3 21 



B. Formal Implementation (Bosonization) 

Formally, the above discussion is equivalent to the 
statement that the low-energy properties of the 1DEG 
are directly related to the properties of two independent 
bosonic field theories with Hamiltonian densities 

= y [(d x 9 a ) 2 + (d x cb a ) 2 ] + V a cos(/3 Q Q ) (1) 

where a — c, s for the charge and spin fields, respectively, 
9 a is the dual field to <j) a , or equivalently d x 9 a is the mo- 
mentum conjugate to <p a . V a can be set equal to zero 
PH at the gapless fixed point. Otherwise, when V a is 
non-zero and relevant, i.e. when f3 a < V&tt, it sets the 
scale of the gap, |2p according to the scaling relation 
A Q ~ V a (A/V a ) 1 - l3 <x /8n where A is an ultraviolet cutoff 
parameter. The Luttinger exponents, K a , determine the 
value of p a to be y/8nK a , and also specify the relation- 
ship between correlation functions expressed in terms of 
the bosonic field operators, and physical correlation func- 
tions, expressed in terms of the original electronic field 
operators, 



1pX,tr( x ) = Mr CX P [iXk F X - i3>\ ja -(x)] 



(2) 
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where M a contains both a normalization factor (which 
depends on the ultraviolet cutoff) and a "Klein" factor 
jlO| (which can be implemented in many ways) so that 
V'a.o-^) anticomutes with J\f a > for a ^ a' and commutes 
with it for a = a 1 . Also 



$A, CT = 



(3) 



where A = ±1 refers to left and right moving electrons, 
and a = ±1 refers to the spin polarization. In Eq. (|3j) , we 
have expressed the fermion operators in terms of "bare" 
bosonic fields, 4> a , which are related to the interaction 
shifted normal fields |l0| that appear in Eq. (|l|) by the 
canonical (Bogoliubov) transformation, 



(4) 



This transformation brings the Hamiltonian into canon- 
ical form, so that the Luttinger exponents appear only 
in the relation between the fermionic and bosonic fields, 
and implicitly in the values of (3 a . 

;From Eq. (§), it is a straightforward (and standard 
fUjl ) exercise to obtain bosonic representations of all in- 
teresting electron bilinear and quartic operators. Physi- 
cally, (j) c and 4> s are, respectively, the phases of the 2kp 
CDW and SDW fluctuations, and 9 C is the supercon- 
ducting phase. The long-wavelength components of the 
charge (p) and spin (S z ) densities are given by 

P (X) = V'l.aV'A.a = 2fc F /7T + y/2K c /n(d x <l> c ) 
S g (x) = (l/2)^(7^ )CT VA,a = ^2K s /n{d x 4> s ). 

A,<T 

We also explicitly bosonize the 77-pairing operator, |l6| ] 
whose correlations are sensitive to the presence or ab- 
sence of a charge gap, 



Vx = ^1,^1,1 



(5) 



exp[iy/2iv/K c (0 c + K c cf> c ) + 2ik F x}. 



This operator is not usually studied as its scaling dimen- 
sion S v = (K c + K ( T 1 )/2 is greater than 1 for K c > 0, 
so the corresponding susceptibility is never divergent. 
However, it is interesting in the present context as it 
has finite momentum and is independent of the spin 
fields. It is easy to see from Eq. (j^) that spatial trans- 
lation by xq is equivalent to the chiral transformation 
<j>c —> <j>c — kFXo t/2/itK c , which must be a symmetry of 
the Hamiltonian for an incommensurate system. Simi- 
larly, gauge invariance implies that the Hamiltonian is 
invariant under 6 C — ► 6 C + const. As a consequence, the 
Hamiltonian must depend only on derivatives of 9 C and 
4> c , so the r\ operator defined above must always create a 
gapless excitation. 

A final comment is in order at this point. The abelian 
representation favored in the present paper is not man- 
ifestly spin-rotationally invariant. This is an advan- 
tage whenever spin-rotational symmetry is broken at the 



Hamiltonian level, as there is no need for special treat- 
ment of symmetry-breaking terms. Spin-rotational in- 
variance implies a specific value of the Luttinger expo- 
nent K s , which may be obtained by comparing the spin- 
spin correlation functions for different spin directions. 
For example, where there is no spin gap, spin-rotation 
invariance can easily be seen to imply K s = 1 at the 
fixed point. (Slow flows as K s approaches 1 also can give 
logarithmic corrections to various correlation functions) . 
So far as we know, in all cases studied to date, the fixed 
point value of K s in a spin-gap phase is K s = 1/4 (or in 
other words, f3 s ~ \/2~n) , at which point the spin correla- 
tions are asymptotically equivalent to those of a dimer- 
ized spin-1/2 Heisenberg model, but it is conceivable that 
other discrete values could occur in other circumstances. 



III. THE MULTICOMPONENT 1DEG 

This section begins with a formal bosonized description 
of the multicomponent 1DEG, continues with a discus- 
sion of the generalized Luttinger theorems for this prob- 
lem, and concludes with a detailed analysis of the specific 
example of the two-component 1DEG. In particular, in 
this latter part, it will be shown how the perturbative 
stability of each potential fixed point Hamiltonian can 
be assessed. 



A. Bosonizing the multicomponent system 

First consider a system of N distinct IDEG's, which 
may be bosonized as in Eq. (0): 



^b,xA x ) = exp [i\k b F x + i$x,<j{ x )] 



(6) 



(b) 

A,CT 



+ \~4>^ 



~0? ] + \^ 



where b = 1 to JV labels the different subsystems, the 
Klein factors anticommute for (a, a) ^ (b,cr'), and 
the bosonic fields satisfy canonical commutation rela- 
tions, 



[^ ) (y),djf ) (x)\ = id a , b 5 a , 5(x - y). 



(7) 



The tilde field variables appearing here are the bare fields, 
unshifted by interactions. 

In the continuum limit, the Hamiltonian density of this 
system consists of three terms 



TL = Ti. c + TL S + TLi 



(8) 



Here 7i c includes all the marginal interactions involving 
the charge degrees of freedom, 

n c = i {(dJZ)w c (dJ c ) + {d s 4%)v B (aJ B )} , (9) 
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where 9^ and <p c are column vectors with N components 

9 C i and <j) C i respectively, and W c and V c are real, sym- 
metric N x N matrices. So that the spectrum is bounded 
below it is necessary and sufficient that the all eigenvalues 
of W c and V c be positive. Ti s is similarly defined for the 
spin degrees of freedom; however, at the spin-rotationally 
invariant gapless fixed point, W s = V s = u s , where 
[u s ] a ,6 = u^S a .b is the spin- velocity matrix. Finally, 
TLint contains the terms nonlinear in the field variables 
(the various cosine interactions), which when relevant 
lead to the opening of gaps in the spectrum, and when 
irrelevant can be ignored. 

For the case in which the nonlinear interactions are 
perturbatively irrelevant, TL is the fixed point Hamilto- 
nian for a system with N gapless charge and ./V gap- 
less spin modes. In this case, as in the single-component 
problem, we perform a Bogoliubov transformation to nor- 
mal coordinates. |2j] A more detailed derivation is given 
in Appendix A. First define the column vector rj. such 
that 

W c V cRi = ul r,. (10) 
where the the normal mode velocities, and 



(11) 



With these definitions, it is straightforward to show that 

w - = Eizaf ( 12 ) 



and 



v c = w c - 1 ^^ 2 . ! fw c - 1 . 



where a and b are N component real vectors, and its zero- 
temperature equal-time correlation function is given by 

-28 



(6(x)&(0)) = [A\x\ 



where A is an ultraviolet cutoff, and the scaling dimen- 
sion 



1 r 

4tt 



a T M- 1 a + b T Mb 



where 



with 



N = WWW 2. 



(18) 

(19) 
(20) 



(13) 



The perturbative stability of the free-boson fixed point 
can be readily analyzed by studying the scaling dimen- 
sion of the various operators which enter into Hint- As 
usual, the stability of the fixed point turns on whether 
there are any physically-allowed vectors a and b that lead 
to a scaling dimension less than 1, which would imply 
that the operator is relevant. If any of these interactions 
is relevant, it is necessary to identify the new fixed point 
to which the system flows, and to study its properties. 
Typically, the effect of a relevant interaction in Ti.i nt is to 
freeze out certain fluctuations (i.e. to gap some modes) 
and at the same time produce a renormalization of the 
matrices V and W. This leads to a new fixed point 
Hamiltonian, whose stability must be reassessed, since 
operators that were irrelevant at the original fixed point 
could become relevant at the new fixed point. This sta- 
bility analysis will be performed more explicitly in the 
two-component example discussed below. 



Then the Hamiltonian may be diagonalized by a canon- 
ical transformation to new fields and their conjugate 
momenta d x &, : 



(14) 



In transformed variables, the Hamiltonian consists of N 
decoupled acoustic normal modes, 



(15) 



The relation between the fermionic fields and the normal 
mode coordinates is easily derived from this expression 
and Eq. ©. 

Finally, correlation functions of the untransformed 
fields can be expressed in terms of the transformed fields 
using equation (|l4|). A typical operator has the form 



6(x) = exp[i (a T 4>(x) + b T 6(x)) 



(16) 



B. Generalization of Luttinger's Theorem 

The generalized Luttinger's theorem jl5j imposes an 
important constraint on the allowed momenta at which 
gapless neutral excitations occur. No matter how com- 
plex the system (e.g. no matter how many bands cross 
the Fermi surface), unless there is an even integer num- 
ber of electrons per unit cell, there must be a zero en- 
ergy excited state with charge and (for the case of zero 
net magnetization) with crystal momentum 2k F = ttut 
where tit is the total electron density, including all bands. 
Thus if a multicomponent system has gapless modes at 
only one crystal momentum, it must be 2k F , and in a 

system with multiple values of 2k F ^ , there must be a set 
of integers to,;, such that 



2k\ 



N 



b=l 



m^k^p 



reciprocal lattice vector. (21) 



(If some modes are charged, there is an obvious further 
constraint on the integers implied by the neutrality of the 
composite mode at 2k* F .) 
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There is a second general constraint governing the ex- 
istence of a gapless charge excitation, which is discussed 
here for the first time. 

This argument generalizes our earlier discussion of the 
77-pairing mode. Consider the generalized 77 operator, 
which creates 2N C right-moving electrons with spin 0: 

N c 

6=1 

= ( ) exp[i\/27r(6*c + 4>c) + i2k* F x] 
where N c is the number of "extended" charge modes and 

6c = ^20i b) /^Nc and 4> c = ^$ b) /y/W e (23) 

6 b 

are the global superconducting phase and the dual CDW 
phase. Global gauge invariance implies that the Hamil- 
tonian is invariant under the transformation 8 C — > 6 C + 
const. Operationally, "extended" charge modes are de- 
fined to be those modes that acquire a non-zero phase 
under a global gauge transformation. Similarly, spa- 
tial translation is equivalent to the phase shift, <p c — > 
4> c + const.. As a consequence of these invariances, Jimt 
must dependent only on derivatives of 8 C and (so long 
as the system is incommensurate) on derivatives of <j) c . 
Thus, the associated modes must be gapless. This im- 
plies that 7]T t i (and of course, ?7t,-i as well) must create 
a gapless, spin charge 2N c e excitation with crystal mo- 
mentum 2k F , and that the ^-correlations must fall like 
a power-law with distance. In nicLiiy crises, — 
and N c = N, the number of "bands" which cross the 
Fermi surface, but we will encounter cases, such as the 
Kondo-Heisenberg model discussed in Section IIIB, in 
which N c < N. 

This proof relies on the field-theoretic representation 
of operators; it is desirable to generalize it to the actual 
lattice system, but we have not yet succeeded in doing 
so. 

Since a gapless, spinless, neutral excitation with mo- 
mentum 2k p always exist on general grounds, it need not 
be listed when classifying phases. 

C. Classification of fixed points 

The essential steps in extending the above analysis are 
to identify the possible fixed points of a multicomponent 
system, and then examine their perturbative stability. As 
for the single-component system, states are identified by 
their discrete broken symmetries and by the "irreducible" 
or minimal set of charge and spin carrying gapless exci- 
tations. The gapless, spinless excitation at 2k* F implied 
by the generalized Luttinger's theorem may be left im- 
plicit. In contrast to the single-component 1DEG, it is 
necessary to specify not only the modes but also their 
momenta (aside from 0), which are no longer completely 



determined by the generalized Luttinger's theorem. For 
instance, it will be seen that it is possible to encounter a 
state with a gapless charge excitation at crystal momen- 
tum 2kp , and a gapless spin mode at crystal momentum 
2k' F ^ 2k F . Such a state will be labelled [2k F : c) [2k' F : s] 
or, leaving the values of the crystal momentum implicit, 
[c] [s] . In the canonical ordering to be adopted here, the 
modes with the larger momentum, 2k F > 2k F , will be 
listed first. An interesting feature of the multicomponent 
system which is obviated by the generalized Luttinger's 
theorem for the single component case is that, when there 
are multiple values of hp, their values can (and generally 
will) shift continuously as a function of interactions. 

Of course, it is always implicit that, if there exists a 
gapless excitation at crystal momentum 2k F , then one 
can make gapless excitations with integer multiples of 
2k F , as well. However, it is clearly not sufficient to spec- 
ify the number of gapless modes, as proposed by Lin, Ba- 
lents, and Fisher, jll]] For example, the state with gapless 
charge and spin modes with the same crystal momen- 
tum, [cs], is thermodynamically distinct from the states 
[c][s] and [s][c] in which they occur at distinct crystal 
momenta. 



D. The two component 1DEG 

To make the discussion more concrete, and in par- 
ticular to illustrate the nature of the stability analysis, 
we now consider the case of a two component 1DEG, 
N = 2. Two independent, decoupled, and generally in- 
equivalent IDEG's are separately described by an appro- 
priate fixed point (free boson) Hamiltonian. Clearly no 
coupling between the two subsystems can be generated 
by any reasonable renormalization group (RG) transfor- 
mation. Thus the fixed points may be specified for each 
subsystem separately. The discussion will be restricted 
to the spin-rotationally invariant case, although this is 
easily generalized. 

The next step is to determine the circumstances in 
which each fixed point is stable with respect to weak in- 
teractions between the two IDEG's. In general, whenever 
a given fixed point is stable for some range of parameters, 
there is no more to say. If the fixed point is unstable, the 
character of the stable fixed point to which the Hamil- 
tonian flows under renormalization must be determined. 
The new fixed point could, in principle, have only gapless 
modes, although, usually some modes that were gapless 
become gapped. 

If the two subsystems are mutually incommensurable, 
{k F ^/k F ^ = irrational) the only potentially perturba- 
tively relevant couplings are those that do not transfer 
momentum between the two systems. The most relevant 
terms are quartic in fermion operators and are of three 
types. The interaction piece of the Hamiltonian density 
is given by 

H' = H[ +H' 2 +H' 3 (24) 
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Here f are the Pauli matrices. Of these interactions, 
the charge and current-density interactions in 7i' 1; are 
marginal, i.e. they are quadratic in boson variables, and 
so must (and can) be absorbed into the definition of the 
fixed-point Hamiltonian density and treated exactly, as in 
subsection A. The perturbative stability analysis is then 
performed with respect to the remaining interactions, 7i' 2 
and H 3 , by computing the scaling dimensions of these 
operators, as in Eq. (|l8|). 

This stability analysis was carried out previously for 
the two chain problem in Refs. and jq|. The results 
are algebraically complicated, but are simplified, without 
significant loss of physical insight, by considering systems 
in which 



V' = -(v^/v^)(kPkP)V, 



(34) 



where kP and vP are the Luttinger exponent and ve- 
locity at the decoupled fixed point of subsystem 6=1, 
2, with a — c, s for charge and spin modes, respectively. 
In this case, 



\(kP + kP + i/kP + i/kP) 



(35) 



^(Ai a 2 ) - \ (a/kP + b/kP + K p + K p] 
<Wa 2) = \ (a/kP + b/kP + i/ K p + i/ K p) , 
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AW 



(36) 



B= 1- 



2VK, 



nvP 



U 1 + 



4VV 



(nvP) 



The weak-coupling fixed point of two incommensurate 
IDEG's is stable if all of these dimensions are greater 
than 1, and is unstable otherwise. For weak interactions, 
any gapless charge modes have kP near 1, and kP 
generally increases with increasingly strong repulsive in- 
teractions. When the fixed point Hamiltonian is spin 
rotationally invariant, K S j = 1. 

The expressions in Eq. (35) were derived for the gap- 
less fixed point Hamiltonian, but it is relatively easy to 
deduce how these scaling dimensions are altered at a 
strong-coupling fixed point in which certain fluctuations 
are frozen out by the presence of relevant interactions 
of the form cos[f3<f)P] and/or cos [f3<j) P] which open a 
gap. This is equivalent to replacing kP by an effective 



Luttinger exponent, K, 



(b) 



and/or K. 



(b) 



0. Con- 



versely, if the fluctuations of the dual phases are sup- 
pressed by a relevant interaction of the form cos[{39P] 
and/or cos[f39P], these expressions should be evaluated 



in the limit kP — > 00 and/or Kf 1 — > 00. Other types of 
strong-coupling fixed points can be analyzed in the same 
fashion. 

It is worth commenting, briefly, on the physical impli- 
cations of the dependence of these various scaling dimen- 
sions on the parameters in the fixed-point Hamiltonian: 

The spin interactions, H.' 2 , are manifestly unimpor- 
tant {$(^.§2) and < ^(a 1 a 2 ) are mnn ite) if either subsys- 
tem has a spin gap. This makes good physical sense. 
If neither system has a spin gap, then the constraints 
of spin-rotation invariance imply that these interactions 
are marginal {^g 1 .§ 2 - ) — 1); further analysis (i.e. carry- 
ing out the perturbative analysis to order J 2 ), follow- 
ing on the work of Sikkema et al. on the Kondo- 
Heisenberg problem, shows that for antiferromagnetic 
couplings (J > 0), these interactions are marginally rele- 
vant while for ferromagnetic couplings (J < 0), they are 
marginally irrelevant. Sikkema et al. speculated that in 
the antiferromagnetic case, the system scales to a strong- 
coupling fixed point with J large and a total spin gap. 
This conclusion is supported by numerical studies carried 
out by these same authors, and by additional analytic 
work by one of us. |l7j 

The singlet pair tunnelling interaction, Tt' 3 , has its scal- 
ing dimension significantly reduced if either or both sub- 
systems have a spin gap, since then the effective kP in 
Eq. (35) is zero. For instance, if subsystem 1 has a spin 
gap, and subsystem 2 does not, then from Eq. (35) 



(b) 



'(AiA 2 ) 



\{a/kP + b/kP 



1 

4' 



(37) 
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The underlying physics is analogous to the proximity ef- 
fect in superconductivity, and we have named it [|| the 
"spin gap proximity effect:" because subsystem 1 has 
a spin gap, it is already substantially superconducting, 
pq| so it can readily infect any coupled subsystem with 
its superconducting character. From this point of view, 
one would expect a relevant pair tunnelling interaction to 
induce pairing correlations in subsystem 2 (i.e. to open 
a spin gap) and to lock the superconducting phases of 
the two subsystems (i.e. to gap the out-of-phase CDW 
mode). We have confirmed the validity of this expecta- 
tion by exact an solution of this problem in a particular 
solvable limit. [|| When J s is relevant, i.e. whenever 
<5(Ai-A 2 ) < 1) the stable fixed-point behavior is charac- 
terized by a total spin gap, and a locking of the charge 
degrees of freedom of the two subsystems. Indeed, this 
effect is very efficient at destabilizing any fixed point with 
a partial spin gap; for instance, any [cs][c] fixed point is 
most often unstable for non-zero J s , and flows to the [c] 
fixed point. The effect of a relevant triplet pair tunnelling 
interaction, J t , has not yet been thoroughly investigated. 

When the two subsystems are mutually commensurate, 
or nearly commensurate, the above stability analysis be- 
comes more complicated. We defer detailed discussion of 
this problem to a later date. However, a few interesting 
features of the problem can be understood on the basis of 
very general considerations. In the first place, the decou- 
pled fixed point of two IDEG's with the same values of 
2Uf is always perturbatively unstable, unless at least one 
system has a fully gapped spectrum. This follows directly 
from the observation that in all the known phases of the 
1DEG, at least one susceptibility is enhanced relative to 
non-interacting electrons. What this also implies is 

that, if we start with two decoupled IDEG's with 2k^p 

(2) 

nearly equal to 2k F , and then gradually turn on interac- 
tions between them, there is a strong tendency to induce 
transfer of electrons between the two subsystems, with a 
cost of unperturbed energy but a gain of energy from the 
relative commensurate locking of the two subsystems. As 
a result, one expects a relative incommensurate to com- 
mensurate transition as a function of increasing interac- 
tion strength in such systems. As for the 1DEG itself, the 
situation is somewhat more complicated for higher order 
relative commensurabilities, since newly-allowed interac- 
tions are generally irrelevant when the interactions in the 
1DEG are weak, and relevant only when they are suffi- 
ciently strong and sufficiently long-ranged. P6| 

IV. APPLICATION TO SPECIFIC MODEL 
PROBLEMS 

To demonstrate the utility of this analysis, we conclude 
with a discussion of four specific problems that have been 
of considerable recent theoretical interest. 



A. Critique of the Perturbative RG Analysis of the 
N-chain Hubbard model 

There have been a number of recent papers concerning 
the phases of the N-chain Hubbard model, following the 
early work of Varma and Zawadowski. |27| In particular, 
in two interesting papers, Balents, Fisher and Lin [Tl| ] 
analyzed the renormalization group flows in the neigh- 
borhood of the non- interacting fixed point, by computing 
the beta function to lowest order in powers of U/t, where 
U is the on-site repulsion between electrons and t is the 
intersite hopping matrix. (More generally, they allowed 
for possibly different values of the hopping amplitudes, 
t and t' ', respectively parallel and perpendicular to the 
chain direction.) Specifically, they identified which in- 
teractions are perturbatively most relevant for various 
geometries of the chains, and as a function of the elec- 
tron concentration per site and the ratio t' /t. They then 
conjectured a phase diagram by analyzing the nature of 
the fixed point obtained by bosonizing the model with 
the relevant interactions taken to infinity, and all others 
neglected. 

This analysis has, we believe, three flaws, which lead 
to significant errors in the resulting phase diagram and 
other conclusions. 1) The only rigorous conclusion that 
can be drawn from a perturbative RG analysis when 
there are relevant interactions is that the initial fixed 
point is unstable, and that therefore the asymptotic 
physics is controlled by another fixed point. (Tracing 
the effects of the perturbatively relevant interactions to 
strong coupling is, at best, suggestive of the character 
of the new fixed point.) 2) Even if we accepted the 
perturbative analysis of the nature of the interactions 
which are important at the strong-coupling fixed point 
to which the Hamiltonian flows, it is essential to perform 
a perturbative stability analysis at the new conjectured 
strong-coupling fixed point to make certain that it is, in 
fact stable. This was not done by Balents, Fisher and 
Lin. Specifically, their analysis ignores the possibility 
that interactions that are marginal or irrelevant at the 
weak-coupling fixed point, can become relevant at the 
strong-coupling fixed point. 3) Because of the presence 
of marginal interactions, e.g. interactions which shift the 
various collective mode velocities, it is an incorrect pro- 
cedure to perturb about the non-interacting fixed point. 
The proper perturbative RG analysis should include the 
effects of the marginal interactions exactly, as in Sec- 
tion IIIA, and should be performed about the appropri- 
ate free boson fixed point. Since the various collective 
mode velocities enter the one-loop RG equations, this 
makes small (but at times important) differences in the 
character of the weak-coupling flows 
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1. Stability of Phases: 

How do these general observations affect the conclu- 
sions of Lin, Balcnts, and Fisher concerning the 
phase diagram? Where they find no relevant interactions 
in their approach, and conclude that all of the modes of 
the non- interacting system remain ungapped {e.g. where 
they find phases of the form [cs] n , or CnSn in their no- 
tation), or where they find completely gapped phases, or 
phases with only a single, ungapped charge mode (e.g. 
phases labeled [c] or CI SO in their notation) the only 
differences between our analysis and theirs are 0(U/t) 
shifts of the locations of various phase boundaries due to 
the effect of marginal interactions. However, all of the 
partially spin-gapped phases, (such as the phase [cs][c], 
which, in their notation is C2S1), are destabilized by the 
spin-gap proximity effect. To see this, note that the pair 
tunnelling interaction, like all other interactions, has scal- 
ing dimension 5 — l+0(U/t) at the non-interacting fixed 
point; the opening of a spin gap, with all other interac- 
tions held small (of order U/t) reduces this dimension by 
a finite amount, e.g. to 5 = 3/4 + 0(U/t) in the exam- 
ple in Eq. (|37|). Thus, of all the conjectured phases in 
their phase diagrams, only the familiar phases with one 
or fewer gapless charge and spin modes, and the novel 
totally gapless phases, [cs] m , are stable at weak interac- 
tions. 

Of course, as is implied by point 1) above, it cannot 
be ruled out that during the flow to strong coupling, in- 
teractions other than the ones identified in perturbation 
theory will get large, and this could alter the stability of 
the partially-gapped phases. However, we suspect that 
the partially-gapped phases are not, generically, stabi- 
lized for large U/t, because strong interactions can pro- 
duce significant shifts in the values of kp\ which then 
permits the commensurate locking of different subsys- 
tems; indeed, this conclusion was reached, previously, by 
Schulz. |2^] In this regard it is worth noting that there 
have been extensive numerical | p9| studies of various N- 
leg Hubbard ladders for intermediate to strong U/t, in 
none of which has evidence for these exotic phases been 
reported. 

B. Example of the spin-gap proximity effect: the 
asymmetric 3-chain Hubbard Model 

In Ref. ( |^1), where the spin-gap proximity effect was 
first elucidated, the simplest model system to which it 
was applied was the asymmetric three-leg Hubbard lad- 
der. The predictions made there were later confirmed in 
numerical experiments on the symmetric three leg t — J 
ladder. (3) For the reasons outlined above, these results 
are in clear disagreement with the predictions based on 
the perturbative RG analysis of Lin, Balents, and Fisher, 
jllj We briefly review the analysis here, as an illustrative 
example. 



As in Ref. ( (|]), consider an asymmetric system, with 
a two-leg Hubbard ladder weakly coupled to a one-leg 
Hubbard "chain" , with a difference in site energy e. For 
concreteness, take U to be large, and consider the phase 
diagram as a function of electron density, although in Ref. 
( ll) it was constructed as a function of e. The method 
of analysis, as presented above, first neglects the coupling 
between the ladder and the chain, and then assesses its 
effect on the final result. 

When the electron density per site is n — 1, the system 
is manifestly insulating. The two-leg ladder, on its own, 
has a spin gap while the chain has a gapless spinon mode. 
Because of the spin gap, this fixed point is stable for weak 
chain- ladder couplings, so the phase is [s]. 

For e > 0, when the density of electrons is reduced 
slightly to n = 1 — x with the number of "doped holes" 
x « 1/3, the added holes go onto the chain. Be- 
cause the interactions in the chain are repulsive and spin- 
rotationally invariant, and the electron density is incom- 
mensurate, the chain will form a Luttinger liquid state 
with gapless charge and spin modes. Pair tunnelling be- 
tween the chain and the ladder induces an effective at- 
traction between spinons, but because there is an energy 
denominator It* (where e* is the renormalized energy to 
transfer a singlet pair of electrons from the ladder to the 
chain) , the bare repulsion between electrons on the chain 
is the dominant interaction, so the decoupled fixed point 
is stable; this phase is [cs] . 

Finally, with increasing doping, although still in the 
regime x < 1/3, the value of e* decreases steadily due 
to the repulsion between doped holes so, if |e| is not too 
large, we reach a regime in which the pair tunnelling 
between the chain and the ladder becomes significant. 
Now, via the spin-gap proximity effect, the chain becomes 
infected with the spin gap of the ladder. The result is a 
phase which has a total spin gap, and only the charge 
2e and neutral spinless modes implied by the generalized 
Luttinger's theorem; this phase is [d. 

The system studied numerically j| is, in fact, the three 
chain t — J ladder. The differences between the t — J 
and Hubbard models for intermediate to strong U are 
not believed to be very significant in the present con- 
text. Because of the boundary conditions, the central 
chain of the ladder is physically distinct from the two 
edge chains. Even though the bare difference in site en- 
ergy e = 0, there is manifestly a non-zero value of e*. 
Thus, there is, in fact, a very close (although not entirely 
quantified) relation between the system we analyzed the- 
oretically, and that studied in the numerical experiments, 
so it is not surprising that the reported phase diagrams 
agree. In the numerical experiments, with J/t ~ 0.5, 
the critical value of doping at which the transition from 
[cs] to [c] occurs is x ~ 0.06; the very small value of this 
critical density reveals the robustness and strength of the 
spin-gap proximity effect. The perturbative RG analysis 
of Lin, Balents, and Fisher leads to a phase diagram in 
which the undoped system is (correctly) in the [s] phase, 
but has the doped system exhibiting the [s] [cs] phase over 







the relevant range of x; this is a specific case of a partially 
spin-gapped phase which, as we argued above, should be 
generally unstable to the formation of a fully spin-gapped 
[c] phase due to pair-tunnelling interactions. 

C. Concerning the N — 6 t — J Cylinder 

Recent important advances in the numerical evalua- 
tion of the ground-state properties of correlated systems 
have allowed the study of much larger t — J and Hubbard 
clusters than before. White and Scalapino, consid- 
ered 6-component t — J systems with cylindrical bound- 
ary conditions, i.e. periodic boundary conditions in the 
finite direction and open boundary conditions along the 
chains. To draw conclusions concerning the 2d t — J 
model, it is necessary to perform a two-dimensional fi- 
nite size scaling analysis of these results in order to ex- 
trapolate to the thermodynamic limit. So far, it has not 
been possible to do so, and the conclusions of White and 
Scalapino disagree with those of other studies of compa- 
rably large systems |3C]] which did do a finite size scaling 
analysis. Consequently, it is still unclear to what ex- 
tent these results are representative of the actual ground 
state of the two-dimensional system. However, we may 
imagine that the results of White and Scalapino are rep- 
resentative of the ground-state properties of an infinite 
length 6-lcg cylinder, which itself is an example of a mul- 
ticomponent one-dimensional system. 

The principal finding of White and Scalapino (ob- 
tained for J/t = 0.35) is that the ground state ex- 
hibits "stripe" correlations in the expectation value of the 
charge-density operator. For the 6-leg cylinder at small 
density of doped holes, x, the period of the observed den- 
sity oscillations is A6 = 2/3x. (Here, units of length are 
chosen so that the lattice constant is equal to 1.) Of 
course, since this is one dimension, the density-wave or- 
der observed by White and Scalapino on finite length 
systems should be interpreted as the period of power- 
law CDW correlations in the infinite system. Now, the 
value of 2k* F — G + 2ir/\* , where G is a reciprocal lat- 
tice vector, corresponds to density-wave correlations with 
wavelength A@ = l/3a;, so the period found by White and 
Scalapino is twice that required by the generalized Lut- 
tinger's theorem. In other words, the fundamental gap- 
less, spinless neutral CDW mode of the system occurs at 
a wave number (1/2)2tt/A*, and the excitation at 2k* F is 
thus a second harmonic. 

It is also easy to see from the present analysis that 
these cylinders are not good candidates for high temper- 
ature superconductors, since there are no gapless charge 
2e excitations of this system. At present, it is not clear to 
us whether the system supports gapless excitations with 
charge 4e, corresponding to the injection of an additional 
"stripe" into the system, or whether because there is a 
tendency for spin correlations to suffer a 7r phase shift 
across a stripe, it is necessary to inject charge 8e corre- 



sponding to a pair of "stripes." 

D. Kondo-Heisenberg array 

The Kondo-Heisenberg model is the simplest exam- 
ple of a metallic system (here, a 1DEG) coupled to an 
insulating antiferromagnet (here, a spin 1/2 Heisenberg 
chain). The large (or infinite) charge gap in the spin 
chain implies that charge-transfer interactions, such as 
pair tunnelling, are unimportant. The dominant interac- 
tions involve the spin density, i.e. they are the Kondo 
interaction between the conduction electron spin s(x) 
and localized "impurity" spins f), 

H K = J K J2T j -s(x j ), (38) 

3 

where Xj 3X6 the positions of the localized spins and the 
Heisenberg interaction between nearest-neighbor local- 
ized spins, 

H«™ = J H J2rj-r j+1 . (39) 

3 

Also it will be assumed that Jh *C Ep, which is typi- 
cally true in physical applications. The resulting Kondo- 
Heiscnbcrg Hamiltonian is 

H = Hl DEG + H^ els + H K (40) 

where, the subscript refers to the Hamiltonian of the 
decoupled system with 

Hl DEG = -iv F A / dx^d x ^ a . (41) 

a. A 

It will be assumed that the relative concentration of lo- 
calized spins is c — 1/6 < 1, i.e. Xj = jb, that the two 
systems are relatively incommensurate, and that 2k F is 
incommensurate with the underlying lattice. The effec- 
tive Fermi wave number (in the sense of the generalized 
Luttinger's theorem) for the 1DEG and the spin chain 
are 2kp and 2kU els — ir/b respectively, so for the cou- 
pled system, fll5j ] there must be a gapless neutral excita- 
tion with wave number 2k* F = 2k F + f , and a charge 2e 
excitation with momentum 2k F ■ 

The determination of the phase diagram of this model 
provides a further example of the application of the meth- 
ods developed above. We shall give a brief physical de- 
scription of the origin of the various phases — the theo- 
retical manipulations may be found in Appendix II and 
in the published and unpublished literature. Pf7|,p7[| 

1. The decoupled Luttinger liquid: Jh —Jk > 

As discussed above, spin-rotation invariance implies 
that, to lowest order, the spin coupling between two gap- 
less systems is marginal, while to second order, as pointed 
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out by Sikkema et al, |31| the interactions are perturba- 
tively irrelevant for Jk < (ferromagnetic interactions) 
and relevant for Jk > 0. Thus, for ferromagnetic Kondo 
coupling, the decoupled fixed point is perturbatively sta- 
ble. This phase has, trivially, one gapless spin excitation 
at momentum 2kp els , and gapless spin and charge excita- 
tions at momentum 2Uf- Since 2fcp > 2kp els , this phase 
is labelled [c, s; 2kp][s; 2kp els ], and has no discrete bro- 
ken symmetries. In short-hand notation, the decoupled 
fixed point is classified as [c, s][s]. However, for future 
reference, it is important to note that there is also a gap- 
less, odd parity charge 2e composite pairing excitation 
at 2kp els , as discussed in Appendix B. Because this is a 
composite excitation, its existence is already implied by 
the existence of the other gapless modes. However, as 
we shall now show, when the Kondo coupling produces a 
spin gap, this composite mode can still remain gapless. 

2. Antiferromagnetic Kondo coupling: Jk > 

a. Odd parity singlet pairing: Jh 3> Jk > 0. If 
Jk > 0, the decoupled fixed point is unstable, and the 
low-energy physics is governed by a strong-coupling fixed 
point with a spin gap. ]3l[ | This phase has several unex- 
pected features. Of course, as required by the generalized 
Luttinger's theorem, there is a neutral, spinless gapless 
excitation of this system with a minimal momentum 2k* F , 
as pointed out by Yamanaka et al. |l5| Remarkably the 
effective Fermi sea knows about the localized electrons 
as well as the itinerant ones. In addition, this system 
clearly has a gapless spinless charge 2e excitation cre- 
ated by the 77-pairing operator of the 1DEG; in this case, 
since the localized spins are unaffected by a global gauge 
transformation, our theorem implies that this mode car- 
ries momentum 2kp . Of course these two statements, 
taken together imply that there exists a gapless, spinless 
charge 2e excitation with momentum 2kp els . Because 
4kp ms — 2-K/b is a reciprocal lattice vector, 2kp els and 
—2kp ets are equivalent; as a consequence, excitations 
with momentum 2kp ezs can simultaneously be charac- 
terized by their parity. As discussed in Appendix B, it 
may be shown that, in the present case, the only gap- 
less charge 2e excitation with momentum 2kp els has odd 
parity. |32[ This phase is labelled [c][c], or more fully 
as [c; 2kp][(c, odd); 2kp els ] (where we put the modes in 
parentheses when an additional descriptive element, such 
as even/odd parity, must be noted). 

There are some remarkable features of the charge fluc- 
tuations in this state. Q Whereas the decoupled system 
had gapless CDW modes at 2k F and 2k p = 2k F + 2kp Iets , 
the spin gap phase retains only the composite CDW mode 
at 2k* F . (See Appendix B.) As a result, because of the 
mismatch between 2kp els and 2k F , there are no spin- 
less charge 2e gapless excitations at momentum zero, 
and the system cannot be a conventional superconductor 
with pairing at total momentum zero induced by spin 



fluctuation exchange. An intuitive feeling for the ori- 
gin of the spin gap can be obtained by considering the 
strong-coupling limit of the model, although care is al- 
ways needed in identifying the specific strong-coupling 
fixed point to which a given weak coupling Hamiltonian 
flows. In the present case, the natural candidate is a 
model in which the conduction electrons form singlets 
with the localized spins, and any remaining localized 
spins form singlets with each other. |n| 

b. Staggered pairing: Jk ^> Jh > 0. An exact so- 
lution of the model also with a spin gap, may be 
obtained from a field theory with explicitly broken spin- 
rotational symmetry and in a special "Toulouse limit," 
in which one component of the Kondo coupling takes a 
specific value ~ Ep. The renormalization group strategy 
behind this solution was described in the introduction. 
The long-distance behavior of the system in the Toulouse 
limit is spin-rotationally invariant, |l3|,|33| which implies 
that spin-rotational symmetry breaking terms are irrele- 
vant at the fixed point, so it is unlikely that the behavior 
we found is an artifact of the model. This phase may 
be distinguished |52| from the weak-coupling (Jh Jk) 
spin-gap state by classifying its gapless excitations. 

Clearly, as before, there must exist a gapless neutral 
excitation with a minimal momentum 2k F , a gapless 
charge 2e excitation, produced by the 77-pairing opera- 
tor, with momentum 2k F , and as a consequence of these 
two general statements, a gapless charge 2e excitation at 
momentum 2kp els . However, we find Q that there exist 
both even and odd parity gapless charge 2e excitations at 
momentum 2kp els , so that there is, in fact, one more fi- 
nite momentum gapless charge mode in this state than in 
the weak coupling spin-gap state. This phase is labelled 
[c][cc], or more fully as [c; 2kp] [(c, odd) (c, even); 2kp els \. 
This solution also provides an example of the fact that 
an analysis of the relevant operators at weak coupling 
does not necessarily tell us the character of the stable 
fixed point. It is long-distance physics (forward scatter- 
ing) that destabilizes the weak-coupling fixed point, but 
the character of the strong-coupling fixed point is deter- 
mined by short-distance physics. 

An elaborate comparison |32| of the two spin-gap 
phases reveals that the additional gapless mode in the 
Toulouse limit phase may be associated with an addi- 
tional hidden broken translation symmetry. A further 
distinction between the states may be made by consider- 
ing the origin of the spin gap and its consequences for en- 
hanced pairing correlations; an intuitive, strong-coupling 
picture of the origin of the spin gap in the Toulouse limit 
phase involves pairing of the spins in each subsystem, sep- 
arately. The existence of the two spin-gap fixed points 
of the one dimensional Kondo-Hciscnbcrg model under- 
scores the need to consider the explicit solutions of the 
strong coupling fixed points, which do not follow from 
simply establishing the existence of a spin gap based 
on the weak coupling perturbative renormalization group 
analysis. 
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c. Strong coupling: Jk 2> Ep 3> Jh- Finally, for 
completeness, it is important to remark that there ex- 
ists a direct "strong-coupling" limit of the model Jk ^> 
Ep S> Jh, which has distinct physics from either of the 
spin-gap phases discussed above. In particular, via a 
mapping to the t — J model, it has been shown that 
in this limit the system is governed by a Luttinger liquid 
fixed point with no spin gap. f|4||35| (Note, if Jh ~ Ep 
there is also the possibility of a strong coupling spin gap 
phase. 0) 
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the critical exponent S is given by 



2 ' fl?) 



(A7) 



from which Eq. (|T^) follows with the aid of Eqs. (|ll| ) 
and (p^), and 



M = w- i j2 u mJ w ~ 1 



(A8) 



This equation may be rearranged to give Eqs. (19) and 
© by using Eq. ©. 



APPENDIX B: GAPLESS MODES OF THE 
DECOUPLED KONDO-HEISENBERG ARRAY 



APPENDIX A: THE MULTICOMPONENT FREE 
BOSON FIXED POINT 



This appendix gives a derivation of Eqs. (|f), ([if), 
@, @, and @. First of all, Eq. @ may be rewrit- 
ten as the eigenvalue equation of a real symmetric ma- 

trix W| V c Wc , with eigenvalues u 2 ci and eigenvectors 

i 

W c 2 r\ .. It follows that 



Wf V c Wf = W e 4 J] u 2 cl Ri vj Wc 



(Al) 



Equation ( |l2|) may now be obtained by premultiplying 

i _ i 

this equation by W c 2 and postmultiplying by W c 2 V" 1 , 

and using Eq. (10). Similarly Eq. ( J13| ) may be ob- 
tained by premultiplying and postmultiplying Eq. (Al) 

by W" 2 ". 

To derive Eq. (|l8|). first use Eq. ( |l2| ) to rewrite Eqs. 
(O) as: 



and 



b T = ufe' ci b T W- 1 r L . 
Then, if Eq. (|l6| ) is written as 



where 



and 



(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



The study of the different stable phases of the Kondo- 
Heisenberg array begins with an analysis of the gapless 
excitations of the decoupled fixed point. From there, as 
usual, we sort the phases by determining which of these 
excitations become gapped, and which remain gapless in 
the presence of the (Kondo) couplings between the 1DEG 
and the Heisenberg chain. In this appendix, then, we 
analyze the system in the absense of any Kondo interac- 
tions. However, since our ultimate goal is to study the 
coupled system, we will consider the character of gap- 
less excitations constructed of composites of the gapless 
modes of the two subsystems, as well as the excitations 
of each, separate subsystem. 

The low energy spin currents of the 1DEG, s (x), can 
be decomposed into two parts; 



s[x) = Ji (a;) + [Hi (x) e i2k " x + H.c] 



(Bl) 



where J\ 



2 Sa,ct : (t' ^X.o^a^' 



4>> 



and 7i i 



\ S a a' V4 a^cr^'^-i,^' are respectively the k = and 
the k — 2kp components of the SDW (charge-0, spin-1) 
of the 1DEG. 

The Heisenberg chain spin current, fj, may be similarly 
decomposed into a k = part, J T , and a finite momentum 
k = ? part, (— l) 3 n T (where ^f- is the reciprocal lattice 



vector of the Heisenberg chain) ; 

fj = J T (xj) + {-l) J n T (Xj) 



(B2) 



As explained in the introduction, we count only the 
number of finite momentum excitations. It follows from 
time reversal symmetry that, for finite momentum, if 
there is a gapless mode at momentum q than there is 
also a gapless mode with momentum — q. We count them 
as one mode. To summarize, the gapless spin-1 excita- 
tions of the 1DEG and the Heisenberg spin chain, and 
the operator whose correlation function is most directly 
sensitive to it are listed in table-B-I. 
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Table — B — I: Gaplcss SDW excitations 



operator 


wave 
number 


m 


2k F 




7T 

h 



(B3) 



The incommensurate 1DEG has one charge-0 spin-0 
CDW excitation with momentum, 2k F , created by the 
operator 



Ocdw 



1 



2 ^ 

A,cr 



\,a- 



(B4) 



The reader may notice an apparent conflict with the 
"prediction" of the generalized Luttinger theorem that 
there must be a gapless CDW mode at 2k* F — 2k F + 
?■. This conflict is resolved by realizing the existence 
of composite-CDW |52],[53) order parameters which are 
formed by combining a spin-1 SDW of the 1DEG with 
a spin-1 SDW of the Heisenberg chain into a composite 
singlet O c -cDW-, 



Oc-cdw — s ■ t 

= Ji ■ J T + Ji ■ n T 

+ \fti ■ J T e t2kFX + H.c 



+ \ni ■ n T e 



i2kpx 



+ h.c] (-iy 



(B5) 



The staggered component, n ■ n T , has momentum 2k* F — 
2kp + f, and is thus the CDW excitation required by 
the generalized Luttinger theorem. There is also a new 
composite-CDW excitation J • n T at k = ?. To sum- 
marize, the non-interacting two-chain system of a Lut- 
tinger liquid and a Heisenberg spin chain has gapless 
CDW modes at three wave vectors (table-B-II) . 

Table — B — II: Gaplcss CDW excitations 



operator 


wave 
number 


Hi ■ n T 


2fc F + f 


Ocdw 


2k F 


Ji ■ n T 


TT 

h 



(B6) 



Note that the composite-CDW excitations at wave vec- 
tors ¥ and 2k p + ? are not independent, since they can 

be related through a multiplication by the 1DEG Ocdw 
(which has wave vector 2k F ). 

The charge-2e singlet pairing modes also require care- 
ful consideration. In addition to the usual k = BCS 
even parity singlet pairing, 



A,T^A, 

corresponding to right and left 



we note also the existence of an ^-pairing mode, ip^ ^tp x , 
at momentum ±2kp 
going singlet pairs. 



As with the CDW modes, in addition to the singlet 
pairing modes of the 1DEG, it is necessary to consider the 
composite singlet pairing, O c _sp, (a product of a triplet 
pairing in the 1DEG with a spin-1 mode of the Heisen- 
berg chain) which turns out to be odd parity |36|]37|], 



O, 



-sp = -« 



(B7) 



where the sum over the spin-indices of the spinors is im- 
plicit. It is decomposed into two momentum components: 
a uniform k = composite odd parity singlet 



and a k = | 



• J T (B8) 
a staggered composite odd parity 



singlet 



df_ a !ff r (x) 



(B9) 



The pairing and CDW modes can be related through 
the 77-pairing modes by using the identities |32] : 



C DW,n even \ - A 
O c -cDW,V° dd — O c -SP- 



where, 



^odd 



4^ E A ^,T^i 



(BIO) 
(Bll) 

(B12) 
(B13) 



The above relations are a manifestation of the fact that 
for each gapless pairing mode there is a corresponding 
gapless CDW mode. Therefore, we adopt the custom 
of dropping the CDW modes from the explicit notation, 
(but their "trivial" existence should be implicitly under- 
stood in any fixed point which has the corresponding 
gapless charge-2e gapless mode). 

In summary, there are two independent gapless, finite 
momentum, charge-2e pairing modes. 



Table — B — III: Gaplcss chargc-2c pairing excitations 



operator 


wave 
number 


V 


2k F 


Oc-SP 


7T 

h 



(B14) 



The gapless modes in tables I-III characterize the non- 
interacting fixed point of a two-chain system consisting 
of a 1DEG Luttinger liquid with Fermi wave number 2kp 
and a Heisenberg spin chain with a reciprocal lattice vec- 
tor j. In our compact notation, which counts only the 
singlet charge-2e and spin-1 charge-0 modes at finite mo- 
mentum, the decoupled fixed point is denoted by 
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[c, s;2k F ][s; 



or [c, s][s] for short. However, as emphasized in the 
above, this description leaves implicit the gapless, odd 
parity charged excitation at momentum n/b. Since this 
excitation is a composite of the excitations already listed, 
it can be omitted in a minimal labelling scheme. But be- 
cause of this, it looks a bit mysterious that there appears 
a a gapless composite pairing mode with momentum n/b 
in the spin-gap states which appear under the influence 
of a relevant perturbation - it looks (incorrectly) as if 
the relevant interaction is generating new gapless exci- 
tations. From this viewpoint, one might be tempted to 
label the decoupled fixed point [c] [c, s] , and to view the 
spin one excitation of the 1DEG as a composite of the 
three other modes, however far that is from the actual 
physics of two decoupled systems. The non-uniqueness of 
the label associated with each state is an intrinsic feature 
of the approach taken in the present paper. 
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